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Exact diagonalization numerical results are presented for a 32-site square cluster, with two holes 
propagating in an antiferromagnetic background described by the t-J model. We characterize the 
wave function of the lowest energy bound state found in this calculation, which has d^2_y2 sym- 
metry. Analytical work is presented, based on a Lang-Firsov-type canonical transformation derived 
quasiparticle Hamiltonian, that accurately agrees with numerically determined values for the elec- 
tron momentum distribution function and the pair correlation function. We interpret this agreement 
as strong support for the validity of this description of the hole quasiparticles. 



PACS: 71.27.+a, 75.10.Jm, 75.40.Mg 

The behaviour of mobile holes in an antiferromagnetic 
(AF) spin background has been a subject of intensive 
study, in part because of its possible connection to high 
temperature superconductivity |^. The simplest micro- 
scopic representation of this problem is the so-called t-J 
model |2||^, which captures the important AF correla- 
tions of an "AF metal" [^. At present only the one hole 
problem is well understood. That is, results obtained 
from a variety of analytical techniques agree with one 
another and lead to the spin-polaron concept for the t- 
J model charge carriers. In particular, analytical work 
based on the self-consistent Born approximation (SCBA) 
Q is found to be in excellent agreement with unbiased, 
exact diagonalization (ED) numerics ^ on a large, clus- 
ter with full square symmetry. 

In this report we focus on two holes in the same t- 
J model. This problem is important because (i) it al- 
lows for the possibility of two-particle bound states, and 
(ii) because it is the simplest problem that allows one 
to study the interactions between charge carriers in an 
AF background. At present minimal analytical work has 
been done on this problem, and only small clusters have 
been treated numerically. So, a complete understanding 
of this problem is far from at hand. 

The aim of our report is twofold. Firstly, we combine 
analytical and numerical results of the two-hole problem 
to provide a comprehensive study of this important sys- 
tem. We show, for the first time, the excellent agreement 
between the results of these approaches. Secondly, we 
argue that the spin-polaron language provides a natural 
description of the physics of the system and a consistent 
way of interpreting of the numerical results. 

Computationally we have managed, for the first time, 
to determine the two-hole ground state for two holes 
doped into the 32-site square cluster used in our previous 
numerical work In the zero crystal momentum sub- 
space we obtain only one bound state, and it has dx2-y2 
symmetry. We have characterized the ground state by 



evaluating a number of important expectation values, no- 
tably the electron momentum distribution function, and 
the (spatial) pair correlation function. 

We have found that the two-hole ground state de- 
rived from an effective quasiparticle Hamiltonian origi- 
nally proposed in Ref. Q may be used to calculate the 
same expectation values that were obtained numerically 
via ED. Further, these quantities are remarkably similar 
to those obtained via ED, giving strong support to the 
appropriateness of this quasiparticle Hamiltonian. From 
the analytical point of view the low-energy physics of 
the two-hole system is generally described as a system 
of moderately interacting spin polarons. This analytical 
work shows that the dominant effective interactions be- 
tween spin polarons come from the short-range nearest- 
neighbour static attraction and spin- wave exchange. The 
purpose of this paper is to present our new ED results, 
and to subsequently provide support for this description 
of the internal structure of the quasiparticles and their 
interactions. 

We consider the t-J model with Hamiltonian 
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where {ij) denotes nearest neighbour sites, cj^, Ci^ are 
the constrained operators. We dope a 32-site square 
cluster with periodic boundary conditions with two holes 
away from half filling. We focus on a realistic value of 
the AF exchange J/t — 0.3, although many other ra- 
tios have been analyzed. The eigenstates with zero total 
crystal momentum are classified according to their sym- 
metry under the square point group C^v We find that 
the d^2_y2 state has the lowest energy, which is consistent 
with calculations using smaller lattices 0. However, the 
binding is rather weak: for the ratio of J/t = 0.3 we find 
that Efj/t = —0.05. For our cluster, the binding energy 
becomes positive in the range 0.2 < J/t < 0.3. 

The spatial distribution of holes may be characterized 
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by the pair correlation function, defined as 
1 
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where is the number of holes, and Nsir) is the num- 
ber of equivalent sites at a distance r from any given site. 
Results are shown in Fig. |l|. 



0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 
0.00 



ED numerical result 
CT analytical result 




1.0 



2.0 



3.0 



4.0 



FIG. 1. The spatial correlation function, C(r), for two 
holes doped into a square lattice described by the t-J model, 
for J/t = 0.3. Our ED results and the analytical results 
(mapped back onto a 32-site square lattice) are shown. The 
solid lines are simply a guide to the eye. 



The electron momentum distribution has also been cal- 
culated: (n„) = (jlql) = (riqi 



where cl^^, 



Cqo- are the Fourier transform of the constrained opera- 
tors, and our results are shown in Table 1. This quantity 
has been studied in detail in smaller systems ||8|; how- 
ever, our 32-site lattice results are important for such 
a study because the important wave vectors along the 
AF Brillouin zone (ABZ) edge [from (0,7r) to (7r,0)] are 
present in our cluster. Notably, in this ground state we 
do not find any evidence of hole pockets at the single- 
hole ground state momentum (7r/2,7r/2). In fact, along 
the ABZ edge, (nq) has a maximum at (7r/2,7r/2). 

The same problem can also be approached analytically. 
Theoretical studies of the t-J model have resulted in a 
clear understanding of the nature of the low-energy ex- 
citations for the system near half filling The charge 
carrier created by a hole introduced in an AF background 
is described as a spin polaron, viz. as a quasiparticle con- 
sisting of the hole and a cloud of spin excitations. The 
wave function of the spin polaron in an AF background 
can be written as 
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in terms of the AF magnon and spinless hole creation 
operators, ajj and /ij^., respectively 

Quite recently a new approach to the two-hole prob- 
lem in the t-J model has been developed Q. It used a 
generalization of the canonical transformation approach 
of the Lang-Firsov type. An effective Hamiltonian for 
the spin polarons, which includes both types of hole-hole 
interactions, has been obtained and its general form is 
given by 



^e// = ^k/ik/ik + i^qa^aq (4) 
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where and ujq are the polaron and magnon energies, 
respectively; Vk,k',q is the "direct" polaron-polaron in- 
teraction; Mk.q is the "bare" hole-magnon vertex; i^k.q 
is the renormalization form factor. 

The two-hole problem for this effective Hamiltonian 
has been considered in Ref. |^ , and the ground state was 
found to be a bound state of d^i^y-i symmetry for all 
J/t > 0.2. In terms of the spin polaron operators, the 
wave function of the d-wave bound state with the total 
momentum P — can be written as 
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The allowed functional form of ensures the d-wave 
symmetry of the ground state. Furthermore, due to the 
fact that the ground state has broken AF symmetry. 



^k+(7r tt) ~ ~^k- Consequently, this imposes that the 
centres of the polarons are always on different sublat- 
tices, which in turn guarantees that the total S"^ = 0. In 
what follows we show that our comparison leads us to the 
conclusion that the large higher harmonics of play an 
important role in determining the behaviour of C (r) and 
(nq). 

The binding energies of the d-wave bound state ob- 
tained numerically by ED and analytically from the 
canonical transformation work, for the experimentally 
realistic J/t = 0.3, are Ej^" /t ~ -0.05 and E^"^ /t ~ 
—0.02, respectively. The reasonable agreement of these 
energies supports the idea that the systems being studied 
represent the same physics. 

Figure |l| shows our results for the hole-hole correlation 
function C(r) in the d-wave bound state. Very similar 
trends are found in both results: they both show that 
about 45% of the time the holes prefer to stay at the 
nearest and next-nearest neighbour distances. That our 
|0) , (3) analytical work produces such behaviour is not incon- 
sistent with our statements regarding the form of A^: 
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the centres of the polarons are indeed restricted to be 
on opposite sublattices, but the holes are almost equally 
distributed on both sublattices, with the maximum prob- 
ability of separation being at \/2. In fact, our analysis 
of the harmonics in shows that about 50% of the po- 
larons in the bound state are separated by one lattice 
constant. Then, the peak of C(r) at r = arises from 
"strings" of length one. The reason for such behaviour 
is in the mentioned hole dressing due to hoppings. Put 
another way, the quasiparticles are something akin to a 
"spin bag" with the moving hole inside the bag. Thus, 
the spin-polaron picture provides a natural explanation 
for the v^-paradox found here and in earlier numerical 
studies 0. 

Table I shows a comparison of the ED and analytical 
electron distribution functions (nq). As we found with 
C{r), both the numerical and analytical quantities show 
similar trends. A full explanation of the ED (riq) be- 
haviour based on the superposition of the effects of (i) 
the d-wave symmetry bound state and (ii) of the inter- 
nal structure of the polarons will be presented elsewhere. 
Here we focus on the features along ABZ boundary which 
are not disguised by kinematic effects and can be di- 
rectly connected with the form of A^. Our analytical and 
numerical results have a maximum at k = (7r/2,7r/2), 
and have minima at (37r/4, 7r/4) and (7r/4, 37r/4). The 
first feature may be understood as an effect arising from 
nothing more than the d-wave symmetry of the bound 
state. The electron momentum distribution function is 
reduced from its half-filled value of 1/2 when holes oc- 
cupy that momentum state. However, the bare hole num- 
ber is proportional to (A^)^, which is identically zero at 
(7r/2,7r/2). Thus, (jiq) must show a maximum at this 
wave vector, so one cannot find hole pockets for a d-wave 
symmetry bound state. 

In future work we will report on the spin correlations 
that are present in this two-hole ground state. Further, 
generalizations of the t-J to forms that are believed to 
better represent the doped Cu02 planes have been inves- 
tigated, and this work will also be discussed elsewhere. 
This work was supported by the RGC of Hong Kong, and 
the NSERC of Canada. 
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TABLE I. The electron momentum distribution function, 
(riq), for J/t — 0.3, obtained from our ED numerics and the 
CT analytical work. 
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(nq) (ED) 


(nq) (CT) 


(0,0) 


0.528 


0.500 


(V4,7r/4) 


0.527 


0.514 


(V2,7r/2) 


0.494 


0.500 


(37r/4,37r/4) 


0.408 


0.429 


(7r,7r) 


0.406 


0.393 


{n/2,0) 


0.521 


0.516 


(Tvfi) 


0.482 


0.495 


{7V,7v/2) 


0.395 


0.451 


(37r/4,7r/4) 


0.466 


0.477 
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